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Abstract 

In this work, we investigate the radiative leptonic decays D7 g s — > y£v {£ = e, jj) 
at tree level within the non-relativistic constituent quark model and the effective La- 
grangian for the heavy flavor decays. We find that the contributions from the three 
Feynman diagrams are all important. With the full calculation, the decay branching 
ratios are of the order of 10 -5 for D~ — ► -y£i> (I = e, \i) and 10 -6 for D~ — > j£u 
(£ = e, /i), respectively. These decays can be measured at B factories and future 
CLEO-C experiments to determine the decay constants fo s and fjj. 
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1 Introduction 

The pure-leptonic decays of heavy mesons are useful to determine the meson decay constants, 
and they are also sensitive to new physics beyond the Standard Model(SM) [1]. But it is 
well known that the decays of into light lepton pairs (see, Fig.l) are helicity suppressed 
by m}/m 2 D : 



Here Vckm is the corresponding Cabibbo-Kobayashi-Maskawa matrix element. In the case 
of D~(D~) decay, it is V cs {V c d)- Fortunately the helicity suppression can be overcome by a 
photon radiated from the charged particles at the cost of the electromagnetic suppression 
with coupling constant a. It is possible that the radiative decays may be comparable or even 
larger than the corresponding pure leptonic decays [2]. 

Several years ago, D. Atwood et al. calculated B ± (Df) — > •yiu in a non-relativistic 
quark model [3], with a large branching ratio. In fact, they only considered one dominant 
diagram, neglecting other diagrams. They found the order of for the branching ratio 
of Df — > •yiu decay. 

Later on, Gregory P. Korchemsky et al. calculated these decays in perturbative QCD 
approach [4]. They gave larger branching ratios. And C. Q. Geng et al. did the calculation 
in the light front quark model [5]. Their branching ratios are rather smaller. 

In this Letter, We will study the radiative leptonic decays D7^ ^tv (£ = e, /i) care- 
fully at tree level, using the non-relativistic constituent quark model, similar to Ref.[3], but 
including all diagrams. In the following section, we will calculate the processes — > ^tv 
(£ = e,/i) in the framework of the constituent quark model (see, for example [6]). In the 
third section, we will compare the result with some of the previous calculations [3, 4, 5]. At 
last we will conclude the calculation briefly. 




(1) 
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2 Model calculations 



We begin with the quark diagram calculation of D~ — > ^tv (£ = e,/i). There are four 
charged particle lines in Fig.l, which correspond to four Feynman diagrams contributing to 
the radiative decays D~ — > ^tv (£ = e,/i) at tree level, as shown in Fig. 2. However when 
the photon line is attached to the internal charged line of W boson such as Fig. 2d, there is 
a suppression factor of m 2 c /my V . Thus we neglect it for simplicity. To be consistent in the 
following calculation, we will always neglect the terms suppressed by the factor m^jm^. 

The decay amplitudes corresponding to the other three diagrams are 



n 



a+b 



= —i\f2G F eV cs c 



Qc i~ 



(Pc ■ P~i\ 



(Ps ■ P 7 ) 



H c = -iV2G F eV cs (^P L s) 



ypi ■ Pj) 



(2) 



As mentioned in the Introduction, we will use the constituent quark model to reduce the 
amplitudes into the 'hadronic level'. In this simple model, both of the quark and anti-quark 
inside the meson move with the same velocity. Thus we have 



Pc = ( m c/ m Ds )p£ , = (m, / m Ds )p fl Da . 



(3) 



We use further the interpolating field technique [7] which relate the hadronic matrix elements 
to the decay constants of the mesons. The decay constant fp for a charged pseudoscalar 
meson is defined by [8]: 

<0\A fl (0)\P( q )>=tf P q fl , (4) 

In the case of D~, we have 

< 0\^ l5 s\D s >=if D y Da , (5) 



The whole decay amplitude for D s — > ^tv (£ = e,/i) is derived from eqs. (2,3,5) by 
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neglecting the terms suppressed by mi/m^ 

(m Ds m Da \ . v a g 
A m s m c J s 1 1 



A = V2eG F V cs 



(e>fp L v). (6) 



In the D s rest frame, the differential decay width [8] is 



Neglecting the mass of light leptons, we get the differential decay width: 



dT aG 2 F \V cs \ 2 I 



2 f2 



dsdi 144vr 2 m 3 Ds (m 2 Ds - s) s 

with 



x s {m 2 Ds — s — i) 2 + x c P , (8) 



^ = ( 3 -^y, Xc =( 3 -2^Y. (9) 

V m s / V m c J 

The s, i are defined as s = (pe + p v ) 2 -, t = (pe + p y ) 2 - Integrating eqn.(8) in phase space, we 
obtain the decay width 

T = ^^ fh ^ [Xs + Xc] - (10) 
Using a = 1/137, m c = 1.5 GeV, m Ds = 1.97 GeV, \V CS \ = 0.974 [8], we get 

T(D S -> 7 £P) = 2.3 x 10~ 17 x ( ^f; t ^ GeV. (11) 

For the lifetime t(D s ) = 0.5 x 10~ 12 s [8], and the decay constant used as fu s = 230MeV 

[3], the branching ratio is found to be 1.8 x 10~ 5 . From eqs.(9,10), we can easily see that 

the decay width is sensitive to the decay constant /|, , and the constituent quark mass m c 

(or m s ). Any changes of the two input parameters, will result in a big change in the decay 

amplitude. Therefore the prediction of branching ratios remain the accuracy at the order of 

magnitude, unless we can precisely determinate the input parameters. 

1 By neglecting terms proportional to m;/m c , we drop the infrared divergence terms, which should be 
canceled by the radiative corrections of the pure leptonic decay D~ — > Iv [9] . 
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It is worthy of considering the differential spectrum for experimental purposes. Deriving 
from eqn.(8), we obtain 

m Ds dT 1 dT 



r dEj r d\~f 



24A 7 (1-2A 7 ), (12) 



with A 7 = E 1 /m,D s - This result is the same as Ref.[3]. We show the photon energy spectrum 
in Fig. 3 as the solid line. This is clearly distinct from the bremsstrahlung photon spectrum. 

The lepton energy distributions are 
1 dT 36 



{x c (l-2A I/ )[2A I/ + (l-2A I/ )lii(l-2A v )] (13) 
1 d\ v x s + x c 

+ x s [2A„(3 - 5A„) + (1 - 2A„)(3 - 2A„) ln(l - 2A„)]} , 

= ^^K(l-2A,)[2A, + (l-2A,)ln(l-2A,)] (14) 
1 d\i x s + x c 

+ x c [2A £ (3 - 5A £ ) + (1 - 2A,)(3 - 2A £ ) ln(l - 2X e )]} , 

where \ u = E u /itid s , A^ = Ee/mD a . We show the neutrino energy spectrum f jf-, and 
the charged lepton (e, fx) energy spectrum in Fig.3 as dashed and dash-dotted lines, 
respectively. Eqs. (13,14) are consistent with Ref.[3], if we consider only the diagram in Fig.2a 
with the photon connecting the strange quark line like the case in that paper. 

The formulas above can be applied to the case of D~ decay i.e. D~ — > ^iv (£ = e, /x), 
directly. We get the decay width easily: 



T(D- - 7 €P) = ^^fl-ml- [x d + x c ] . (15) 



with 



x d = ( 3 ) , x c = ( 3 - 2- 



m d j \ m c 

Using m d = 0.37 GeV, m D - = 1.87 GeV, \V cd \ = 0.22 [8], we get 



r(/r -> 7 £z7) = 2.9 x io~ 18 x ( ] GeV / 16 n 

v ' J \230MeV J K ' 

For t(D ) = 1.05 x 10~ 12 s [8], the decay branching ratio is 4.6x 10~ 6 with the decay constant 
f D - = 230MeV [5]. Again, without the precise determination of the constituent quark mass 
m d , the decay branching ratio is only meaningful at the order of magnitude. 
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In the case of D~ decay, the formulas for the differential spectra of photon and lepton en- 
ergy distribution are the same as the D s decay, except replacing x s with in eqs.(12,13,14). 
Their numerical results are shown in Fig. 4 as solid, dashed and dash-dotted lines for 7, u, £ 
respectively. From Fig. 3 and Fig. 4, we can see that the differential spectra of Df and D ± 
radiative leptonic decays are very similar. Only the endpoints of leptonic energy spectra are 
different. 



3 Comparison with other calculations 

In Ref.[3], D. Atwood et. al. made the calculation within the non-relativistic constituent 
quark model like us. As stated in the introduction part that they just considered the con- 
tribution of the emission of photon from the strange quark, i.e. Fig. 2a, for they made an 
analogy with B~ decay directly. After our careful calculation, we conclude that the con- 
tribution of the Feynman diagram where the photon is emitted from the initial light quark 
is dominant enough to neglect the other diagrams in the case of B~ decay, but not for the 
case of D~ or D~ . It can be seen at Table. 1 that the contribution of the other two diagrams 
corresponding to Fig. 2b and 2c must be considered because the interference among the three 
Feynman diagrams is large and destructive. That is the reason why their branching ratio of 
Br(D~ — > ^y£p)(£ = e,/z) decay [3] is about four times of ours. 

Gregory P. Korchemsky et al. used the perturbative QCD method to calculate B and D 
meson radiative decays. Their result for D~ decay is very large. In fact, the perturbative 
QCD approach [10] is good for the B meson decays since the energy release is very large there, 
but may not be good for the lighter D meson decays. In addition, out result is consistent 
with that of [5] within the light front quark model in the case of D~ — > jeu. It is instructive 
to calculate with various models, and the accuracy of various models will be tested in future 
experiments. 
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Table 1: Decay width with different diagrams and their relative size. Considering only the 
diagram where the photon is emitted from the initial light quark (Fig. 2a), heavy quark 
(Fig. 2b) or lepton (Fig. 2c), we get r a , Tb, T c , respectively. And considering the three 
diagrams together, we get T a+ b +c . 





r a 


r b 


r c 


Ta+6+c 




B~ 


1.7xl(T 18 


5.7xl(r 22 


5xl0- 20 


1.2xl0~ 18 


1.40 : 0.0005 : 0.04 : 1 


D s 


3.4xl(T 16 


8xl0~ 17 


4xl(T 16 


2.3xl0~ 17 


14.72 : 3.47 : 17.32 : 1 


D- 


2.1xl(T 17 


2.7xl0~ 18 


1.8xl0~ 17 


2.9xl(T 18 


7.30 : 0.94 : 6.03 : 1 



4 Summary 



We have calculated — > / ~f£u (£ = e, fi) decay in non-relativistic constituent quark model. 
We included all the three Feynman diagrams, and found that none of them is small in 
decays. We obtained the decay branching ratios of Dj — > ^iv , D~ — > ^y£u(£ = e, /i) are of 
order 10~ 5 and 10~ 6 respectively. Such a branching ratio for the radiative leptonic decays 
can be measured in the two B factories and the future CLEO-C Experiments. 

Eqs. (10,15) indicate that the decay rate of — > ^iv (£ = e, /i) is proportional to /|> w , 
so one can use it to determine the decay constant /d (s) - On the other hand, it is seen that 
these processes can also be used to test the \V CS \ and \V c d\ if /d s and fo are known. 
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Figure 1: Feynman diagram in standard model for D s — > I v decay. 
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Figure 3: Normalized energy spectra of the decay D~ — > 7 iv. The solid line is for the photon 
energy spectrum, the dashed line is for the neutrino energy and the dash-dotted line is for 
the lepton energy spectrum, respectively. 
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Figure 4: Normalized energy spectra of the decay D~ — > ^Iv. The solid line is for the photon 
energy spectrum, the dashed line is for the neutrino energy and the dash-dotted line is for 
the lepton energy spectrum, respectively. 
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